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Abstract 



A Boltzmann equation is given for the early stages of evolution of the 
gluon system produced in a head-on heavy ion collision. The collision term 
is taken from gluon-gluon scattering in the one-gluon approximation. < p± > 
and < > are evaluated as a function of time using initial conditions taken 
from the McLerran-Venugopalan model. 

1 Introduction 

In an earlier paper we studied the early stages of the gluon system pro- 
duced in a heavy ion collision. In that paper we took the initial conditions 
of the gluon system from the McLerran-Venugopalan model assuming that 
all gluons at or below the saturation momentum in the light-cone wavefunc- 
tion of the McLerran-Venugopalan model are freed in the heavy ion collision 
while all gluons beyond the saturation momentum are not freed|0, ^, |^, ^, ^]. 
We then wrote an equation for the rate of change of the typical transverse 
momentum of a gluon in the system as a function of time. This equation 
takes into account, in an average way, the transfer of transverse momentum 
into longitudinal momentum due to elastic gluon-gluon scattering in the one 
gluon exchange approximation. We were able to follow the decrease of the 
typical transverse momentum of a gluon with time up to times on the order 
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of -^exp[canst / ^/a ] at which time the typical transverse and longitudinal 
momenta become equal, and our procedure broke down. We argued that 
this was likely the time at which equilibrium was setting in and our simple 
approach was not capable of following the gluon system into equilibrium. 

In this note we reformulate our previous discussion in a much more general 
way in terms of the Boltzman equation with the collision term taken from 
elastic gluon-gluon scattering in the one gluon exchange approximation. In 
the logarithmic approximation for integrating over the angle of scattering, 
the same approximation used in Ref.l, the Boltzmann equation reduces to a 
Fokker-Planck equation]^, given in (15) and (16), for the number of gluons, 
f, per unit of phase space. The infrared cutoff on very small angle scattering 
is taken from Ref.l although the exact form of this cutoff is not crucial for 
the existence of an equation of the type given in (15). 

There are two explicit parameters appearing in (15). The first of these 
parameters, ?7, which is time-independent is given by the initial conditions 
for the evolution. In our explicit evaluations we take t] from the McLerran- 
Venugopalan model in terms of the saturation momentum Qg. The second 
parameter, given in (13) should be determined self-consistently from 
the solution to (15). In our explicit evaluations we take ?7_i to be given by 
the McLerran-Venugopalan model at very early times (28), while the time 
dependence of ?7_i can be estimated using (32). We have also assumed to, 
in (10), to be time-independent. In Ref.l, we found the time-dependence of 
to to have a negligible effect at early times and we expect the same to be 
true here. However, before attempting to solve (15) more completely, say by 
numerical means, at non early times one must give a specification of to-, in 
(10), analogous to what was done in Ref.l, but now in the context of the 
Boltzmann equation. 

2 The Boltzmann equation 

A derivation of the Boltzmann equation can be found in Ref.7, and can be 
written as 

% + ^-^^f = c (1) 

where C is the collision term. / is the number of gluons per unit of phase 



2 



space after the valence quarks have separated in a high energy, head-on, 
heavy ion colhsion. Thus / = f{p,x,t) and 

N{x,t)^j£pf (2) 

is the gluon density per unit volume. / is a scalar function. We shall assume, 
for simplicity, that / only depends on z, the coordinate corresponding to the 
axis along which the heavy ions collide, but not on the transverse coordinates 
X. Further we suppose that the physics of the heavy ion collision is boost 
invariant, at least for final state particles not having too large a rapidity in 
the center of mass system. If we parameterize the momentum and position 
of a final state gluon by 

= P± {cosh y, V, sink y) = (po, P, Pz) (3) 

Xfjt — T{cosh 77, 0, sink rj) — {t, x, z) (4) 

then / = fip±, T,y — rj) expresses the boost invariance of the distribution of 
particles in the final state. In terms oi T,y,r] 



(| , , . ^ (| . ^ ^^(| - 1^)/ (5) 

or, taking rj — z — 

where we have expressed, at z — 0, f as a function of t,p± and Pz, variables 
convenient for the discussion to follow. 

Now turn to the collision term in (1). It is straightforward to show that 

C^d,jSp2[{d,-d^,)fh]B,, (7) 

where we use the notation = ^, = ^J-, / = /(i,p_L,P2), /2 = f{t,P2±:P2z)- 
is given by 

Bij = 2na'^{ ^^ ^ )L[(1 - v ■ V2)dij + ViV2j + V2iVj)] (8) 



3 



with V = and V2 = Except for the factor (j^^) the Bij in (8) is 
identical to that given in Ref.7 for QED. The factor L represents a (pertur- 
batively divergent) integral over the (small) angle of scattering which we cut 
off as in Ref.l at a minimum scattering angle Om given by 

O'm-J (9) 

with to slowly varying with time for t not too large. We shall treat to as a 
constant in this paper. Thus, 

Using (8) and (10) in (7) along with (1) and (6) leads to 

C = XN^ vif{{t,P) + 2A7V_ie Vp-ivf{t,p)) (11) 

where 

N{t) J d'pfit,p) (12) 
is the gluon number density while 

N,{t) = Jd'pf{t,p)\p\' (13) 

and where 

Defining tN — r],tN_i = r]_i and using (6) and (11) in (1) finally leads to 
an equation of the Fokker-Planck type 

(^ - ^P.)f = A77e Vj/ + 2A77-ie Vp(vf) (15) 

with 

tf{t,p) = f{t,p}. (16) 



4 



3 Properties of the solution of the Boltzmann 
equation 

It appears difficult to solve (15) even with simple initial conditions. Nev- 
ertheless, one can discern general characteristics of the solutions, for early 
times, even without being able to completely solve the equation. Taking the 
integral with respect to d^p on both sides of (15) and recalling that / fd^p — rj 
one finds 

1" = " (1^) 

which means that the number density of gluons N(t) decreases like 1/t. 
Multiplying both sides of (15) by p\ and integrating over all p leads to 




(18) 



— <pI>+2<pI>^ 2Xri^ - 
In (18) and (19) we have used 




(19) 




(20) 



and 




(21) 



In case / does not decrease faster than pj^^, which is the behavior in pertur- 
bative QCD, it may be useful to replace (18) by 



|:<Px>=A^fe- '"-"''P"^''' ) (22) 



where 




(23) 
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and 

Vipjp) = / rfW-- (24) 
J p 

At early times the values of pz are very small at ;z = so that we may set 
fj^i ~ rj-i and ri{p±/p) ~ ?] so that (22) becomes 

^<p±>^-X^V-i- (25) 
In the McLerran-Venugopalan model used in Ref.l 

iV^ - 1 

fs{t,p) = 6{p^)Q{Ql -Pi)c^^ (26) 
immediately after the ion-ion collision leading to 



While at early times 



so that (25) becomes 



37r 

Qs is the saturation momentum. Defining 



TTT < P± >= -c—^Qsi (29) 



~ 1 f 2 

Q = - £p p^f^ = -Q^ 

T] J 3 



one finds 



< pjQ >= 1 - c^e' (30) 

a behavior having an identical form, for ^ not too large, to that found in 
Ref.l. 
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Now turning to (19) it is clear that at early times the first term on the 
left-hand side of that equation and the second term on the right-hand side 
are small compared to the remaining two terms at early times so 

< {V.IQ? >= (31) 

Over what range of ^ can (30) and (31) be used? Since we are working 
in a logarithmic approximation we clearly need >> 1 in order to take 
the collision term to be given by (11). In the McLerran-Venugopalan model 
^ ^ £n 1/a even when t is as small as 1/Qs so the assumption of large ^ 
should be (parametrically) good for all times after the gluons are freed in 
the heavy ion collision. In obtaining (25) from (22) we have used r]^^ ri_i 
and rj{p^/p) ^ rj. These assumptions are good so long as typical values of 
are much less than typical values of p±. We shall come back in a moment to 
determine when \pz\ ~ p±. In solving (25) we have taken ?7_i to be constant, 
and given by (28) in the McLerran-Venugopalan model. We can estimate the 
corrections coming from the time variation of ?7_i by taking 

V-iit)=V-i^j^- (32) 
If we use (32) in (25) the result 

< PJQ >'= 1 - (33) 

emerges. Eqs.33 and 30 are identical when << 1, however, we can expect 
(33) to be reasonable even when is not small. Thus, we suggest that (31) 
and (33) should be good estimates for the typical transverse and longitudinal 
momenta so long as < > remains less than < >^ . (The precise factor 
in front of the ^^^^ in (33) is not expected to be rehable, however, because 
(32) should be expected to be a reasonable estimate though not a precise 
evaluation.) 

From (31) and (33) we see that < p1> and < p^ >^ become comparable 
when ^ = ^1 with given by 
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When ^ > ^1 we can no longer expect (31) and (33) to be indicative of the 
evolution of the gluon system after the heavy ion collision. Eqs.(18), (19) 
and (22) should still be correct as they follow from (15), but when ^ > we 
are no longer able to estimate the terms appearing on the right-hand sides of 
(18), (19) and (22). While we are unable to follow the system, analytically, 
beyond ^ = we presume that typical values of p±_ and Pz remains equal, at 
z = 0, and that the gluonic system continues to evolve toward equihbrium. 
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